function y = func_q_from_p(p, CONSTS)

    eps_a = CONSTS.eps_a;
    %y = -1i*p; 
    y = sqrt_from_negative_num(eps_a-p.^2); %(1-p.^2).^0.5;
    
end